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Abstract 

The integrable XXX spin s quantum chain and the alternating s^, {s^ — 
= |) chain with boundaries are considered. The scattering of their excitations 
with the boundaries via the Bethe ansatz method is studied, and the exact 
boundary S matrices are computed in the limit s,s^'^ oo. Moreover, the 
connection of these models with the SU{2) Principal Chiral, WZW and the 
RSOS models is discussed. 
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1 Introduction and Review 



Quantum spin chains are ID quantum mechanical models of microscopic degrees of 
freedom, namely spins, with the Heisenberg {XXX) model solved by Bethe 0], 0, 
0] being the prototype. Two types of quantum spin chains exist, known as "closed" 
(periodic boundary conditions) and "open". Closed spin chains have been very well 
studied (see e.g. [0-0), whereas open chains are less widely studied, although of great 
interest. A very interesting aspect of these models is that in the continuum limit, they 
correspond to 1+1-dimensional integrable quantum field theories [^]. For example the 
S matrix that describes the anisotropic Heisenberg {XX Z) model in a certain regime 
also describes a massive integrable quantum field theory, namely the Thirring model or 
the sine-Gordon model. Also, the S matrix that describes the bulk scattering for the 
spin s XXX chain coincides with the one of the Wess-Zumino-Witten iWZW) model 
0. Finally, open spin chains and the corresponding 2D field theories with boundaries 
display a rich variety of boundary phenomena which — for integrable models — can 
be investigated exactly. 

During the last years, there has been an increased research interest on integrable 
models with boundaries, especially after the prototype work of Cherednik and Sklyanin 
0, P]. In particular, Cherednik introduced the reflection equation in order to obtain 
boundary scattering matrices for theories in half line. In the presence of boundaries 
in addition to the bulk scattering, the scattering of the particles with the boundaries 
— described by the boundary S matrix — should also be considered. The boundary 
S matrix satisfies a collection of algebraic constraints, namely the boundary Yang- 
Baxter (reflection) equation [0. Sklyanin considered this equation in the spin chain 
framework and he generalized the Quantum Inverse Scattering Method for integrable 
models with boundaries. Moreover, Fring and Koberle [^] obtained solutions of the 
reflection equation for the affine Toda field theories with boundaries, whereas Ghoshal 
and Zamolodchikov [jlO| found solutions of the reflection equation for the sine-Gordon 
model on a half line. Finally, De Vega and Gonzalez-Ruiz made similar calculations 
for the XXX {SU{2)), XX Z (A^^'), XY Z and any SU{n), A^^li open spin chain fll 



In this study the integrable XXX spin s quantum chain and the alternating s^, 
chain with boundaries — obtained by fusion — are explored. As a warm up, both 

models with periodic boundaries conditions are reviewed and their relation with 2D 

quantum field theories is discussed. 

To describe the models it is necessary to introduce the basic constructing element. 



namely the R matrix, which is a solution of the Yang-Baxter equation |T^, 

-Rl2(Al — A2) i?l3(Ai) i?23(A2) = -^23(^2) -Rl3(Al) -Rl2(Al — A2) . (1.1) 

We focus on the special case where the R matrix is related to the spin s representation 
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of SU{2), obtained by fusion |Ti 



■Ok 



w+S^ Wo - W3SI 



where wq = X + ^, = w± = i and S^, S"^ are the SU{2) generators in the 
s representation and act, in general, on a 2s + 1 dimensional space V = C^*"'"^. 
generators satisfy the following commutation relations 
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and for e.g. s = 1 they become. 
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For s = we obtain the well known XXX R matrix. 
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1.1 The XXX spin s quantum chain 

The XXX (and XX Z) spin s integrable chain has been intensively studied by several 
authors in the bulk (see e.g.|jl^-|2l|). To construct the model we derive the transfer 



matrix considering the R matrix related to the spin s representation of SU{2) ( |1.2| ). 
We introduce a mass scale in our system therefore, we derive the transfer matrix with 
inhomogeneities fl, 

t(A) =troTo(A) (1.6) 

where 

^o(A) = i?^^^(A - n)Rllj,_,iX + fi) ■ ■ ■ RlliX - mlli>^ + ^) , (1-7) 

i?ol is given by ( |1.2| ) and acts on Vq V^. The auxiliary space Vq is a two dimensional 
space whereas the quantum space Vi, i = 1, . . . ,2N, is a 2s + 1 dimensional space. 



After we diagonalize the transfer matrix we find the following eigenvalues [14 



A^iX) = ix-n + ts+'-nx + n + ts + ;-fl[ 

^ ^ j=i - + 2) 

+ (A-r]-zs + lf(A + l^-zs + ^)^n (1-^ 



M 



(A - A, 
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The corresponding Bethe ansatz equations are, 



M 



e,{K - ^)''e,{K + ^^)^ = - n e2(A„ - A^) , (1.9) 



where q = 2s and 

A + f 



e„(A) 



2_ 

\ in ' 



The BAE ( |1.9| ) are similar to the ones found in for the lattice analogue of the 
SU{2) PCM. 

In general, if we fuse the auxiliary space as well, we find that the eigenvalues of the 
transfer matrix are 1113, 



Ar(A) = Y.^kix-nfakix + n)'' 

k=0 

n (A-A,-t) (A-A,+z/ + |) 

/i (A - A,- + t{k - 1) + i) (A - A, + + i) ' ^ • ' 



where 



i-i ■ k-i ■ 

afc(A) = n - + 2 + n + + + ^^-^^^ 

m=k ra=0 

Let us consider the special case where I = 2s and A = ±fi + Aq = ±fi — |(2s — 1). Then 
i?^''(Ao) becomes the permutation operator, therefore we can derive a local Hamiltonian 
for the system (up to an additive constant) 

The corresponding eigenavalues follow from ( |1.10| ) and ( 1.12|) 

1 2 Af ^ 

^ ^ "2^S§ (A, + + f )(A, + i-rn - f ) ^^-^^^ 

also the momentum and spin are given by 

^^^ih ^A, + (-)"fi-f) ^ ' 

S' = qN-M. (1.15) 

The ground state and the low lying excitations of the model can be studied, in the 
thermodynamic limit N —>■ 00. In this limit, the solutions of the Bethe ansatz equations 
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are given by the so called string hypothesis, i.e. the solutions of (|1.9| ) are grouped into 
strings of length n with the same real part and equidistant imaginary parts 



A^ + -(n+l-2j), 



j = 1,2, ...,n 



:i.l6) 



where is real. It is known [n]-[18|, that the ground state, i.e. the state with 
the least energy, is the filled Dirac sea with strings of length q = 2s. The low lying 
excitations are holes in the q sea and also strings of length n ^ q. One can study 
the scattering among the low lying excitations of the model and show |T^-|21| that 
the S matrix, as s — oo, coincides with the one of the SU{2) PCM i.e., it is the 
Ssu(2) ^ Ssui2) S matrix [^. The SU{2) scattering amplitudes (for the singlet 
triplet and singlet respectively) are 

S'o(A) = expj- / e~ 

^ J — oo 



duj 



2cosh(f) 



J' 



5+(A) 



X + i 



X-i 



:So{X). 



[1.17) 



In the scaling limit, X << Q, the excitations become massive relativistic particles flS 
with energy and momentum. 



6(A) 



2e^''^ coshTrA, 



p{X) 



2e-^^ sinh vrA 



:i.l8) 



where the momentum p is defined modir for even number of excitations. 



Faddeev and Reshetikhin showed in ||18|] that for finite s, inconsistencies in the 
counting of the states exist. Therefore, the interpretation of the S matrix as Ssg{s) Ssui2) 



somehow fails; Sscis) is the sine-Gordon S matrix p4|, with the triplet amplitude be- 
ing 



S'q{X) = exp 



{- 



sinh(^(2^ 



l)f 



du 



2cosh(|)sinh(2s^ 



UJ 



;i.i9) 



and obviously Ssg{s oo) Ssu{2)- However, Reshetikhin conjectured in that the 
correct S matrix for the spin s chain is the Srsos{s) (g) Ssu{2) matrix which coincides 
with the S matrix of the WZW model at level k = 2s {WZWk) [g2[, §, H- The 
Srsos{s) is the scattering matrix of the RSOS model, and the spin s is related to the 
restriction parameter r of the RSOS model i.e. r = 2s+2. (for a more detailed analysis 
|]). More specifically, the spin s 



see 



H, i 



1 chain has a hidden super-symmetry, 
which is described by the RSOS part of the S matrix 0. 

A key observation is that, as A; = 2s ^ oo, the Sjisosi^) Ssui2) matrix degener- 
ates to a tensor product of two rational matrices, and it coincides with the S matrix of 
the SU{2) PCM without topological term, found in |jT6|, [jl8|, The later comment 
reflects the fact that the perturbed WZWk as — oo, reduces to the PCM without 
topological term {Snsosis oo) reduces to Ssu{2)), as described in ||26[ . 
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1.2 The alternating s^, spin chain 



Alternating spin chains have been originally introduced by de Vega and Woyanorovich 
in ||2^ and they have been also studied in the bulk by several authors (see e.g. 
). We define the transfer matrix of the chain with inhomogeneities, 



t = troTo{X) (1.20) 

where 

To(A) = Rl,j,{\ - fi)i?^27v-i(A + f^) ■ ■ ■ i?02(A - ^^X(A + fi) , (1.21) 

and K is related to the spin (i = 1, 2) representation (|1.2|) . The eigenavalues of the 
transfer matrix, after we fuse the auxiliary space as well, are given by (see also |27]) 



TT (A-A,-t) (A-A, + ./ + f) 

/i (A - A,- + i{k - 1) + i) (A - A, + ifc + f ) ' ^ • ^ 

where 



l-l ■ k-l 



«i (A) = Yl - is' + t: + im) Y[{X + is' +in + -). (1.23) 

m=k ^ n=0 ^ 



The corresponding Bethe ansatz equations have the form. 



M 



e,i(A„ - n)%2(A„ + (])^ = - n e2(A„ - A;3), (1.24) 

/3=1 



where q' = 2s' and — = 1. The BEA (|1.24|) coincide with the ones found by 



Polyakov and Wiegmann for the lattice analogue of the SU{2) PCM with WZW term 



Again for / = 2s' (j = 1, 2) and X = ±Q + Xi = ±Q- ^{2s' - 1) the R'{Xi) matrix 
becomes the permutation operator and hence we obtain a local Hamiltonian, 

and the corresponding eigenvalues are 

^ 2 M 

^ = "2^ 5g (A, + (-)"fi + f )(A, + i-rn - ¥) • ^^-^^^ 

Moreover, the momentum and spin of a Bethe state are given by 
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= (gi + q^)N - M. 



(1.28) 



It has been proved that this model has two types of low lying excitations, 



i.e. holes. The scattering among them was studied, the corresponding 5* matrix was 
computed, for — > oo, and it was shown to be the Ssui2) (8> Ssu{2) plus a non trivial 
left right scattering ^'^.^^(A) = tanh|(A — |). The Ssu(2) ^Ssui2), Slr{X) is also the 
massless 5* matrix of the SU{2) principal chiral model with WZW term at level 1 
{PC Ml), conjectured by Zamolodchikovs It has been also proved [^, that in the 



scaling limit A << ^2, both excitations obey a massless relativistic dispersion relation, 
namely 

e\X) = p\X) ~ e-"^e"^ , e\X) = -p\X) ~ e-^^e'"^ . (1.29) 

These are the energy and momentum of the "right" and "left" movers respectively (see 
e.g. p2[) and the factor e"'^^ provides a mass scale for the system. 



2 Spin chains with boundaries 

After the brief review on the bulk XXX spin s and alternating chains, we are ready to 
study these models in the presence of boundaries. For both models we will derive the 
Bethe ansatz equations and we will compute the exact reflection matrices. To construct 
the spin chain with boundaries in addition to the R matrix another constructing ele- 
ment, the K matrix, is needed. The K matrix is a solution of the reflection (boundary 
Yang-Baxter) equation 0, 

i?i2(Ai - As) Ki(Ai) i?2i(Ai + A2) K2(A2) 
= i^2(A2) RuiXi + A2) Ki(Ai) i?2i(Ai - A2) . (2.1) 

In what follows we are going to use Sklyanin's formalism in order to construct both 
models with boundaries. 

2.1 The XXX spin s open chain 
2.1.1 The Bethe ansatz equations 

Let as first consider the XXX spin s = | quantum chain with boundaries. We define 
the transfer matrix for the open chain ^ 

t{X) = tro K+{X) To{X) K,{X) fo(A) , (2.2) 

To (A) is the monodromy matrix (|1.7D and 



To(A) = Rfoi^ - fi)i?^^(A + n) ■ ■ ■i?^^_io(A - ^)RrNoi>^ + ^) . (2.3) 
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The K (A) = K(A,^ ) matrix satisfies the reflection equation (|2.1|) and the K^{\) 
K-{-X - pY ^ -^+) satisfies, 



RuhXi + A2) K+{X,y^ i?2i(-Ai - A2 - 2p) K+iX^Y' 
KtiX^Y' Rui-Xi - A2 - 2p) K+{X,Y' ^2i(-Ai + A2) , 



(2.4) 



where p = i. We choose the diagonal solution of the reflection equation |]TD|], |TT 
namely 

K{X) = diag{-X + i^, X + i^). (2.5) 
Note that in the presence of boundaries we have to fuse the K matrix as well, and we 



use the following fusion hierarhy for the transfer matrix [^, ||35 

t^{X) = C2,_i(2A+(2j-l)2)[t^-^(A)a(A + (2j-l)z) 
A(A + (2j - 2)2)C2,_2(2A + {2j - 2)^) 



C(2A + 22(2j-l)) 



-f-\X) 



where 



A(A) = A[K+(A)] 6[T{X)] A[K-{X)] 6[f{X)] , 
and the quantum determinants are flj], [Q, 

6[T{X)] = 6[f{x)] = cix + n + ifa^ -n + if, 

A[K-iX)] = g{2X + p)/(A, T), A[K+(A)] = gi-2X - 3p)/(A, t 



(2.6) 



(2.7) 



(2.^ 



moreover 



C(A) 



and 



(A + z)(-A + i), C(A) = {X + is+ '-){-X + IS+'-), 
-A + z, /(A,0 = (A + 2 + 0(A + ^-0 



0(A) = nC(A + ^fc), Co(A) = l. 
fc=i 



(2.9) 



(2.10) 



Then the eigenvalues of the transfer matrix, after we fuse the auxiliary space, are given 
by 

Af^(A) = Y.hkWfkWakix-nf'^a.ix + nY'^ 



(A-A,-f) (A-A, + z; + f) 



fc=0 
M 

n _ ^ _ 1) + I) (A _ A,. + zfc + I) 

il (A + Xj - i) (A + A, + z/ + I) 
/ii (A + A,- + i{k - 1) + i) (A + A,- + iA; + I) 



(2.11) 



where have been aheady defined in ( 1 . 1 1| ) and 

fk{X) = Jt{X)f,{X) 



(2.12) 



where 

fkW = n (-^ + + ^m + ^) Ui-X + t^^ - m) (2.13) 

m=k n=0 

/^(A) correspond to the left (+) and right (-) boundary. hk{X) are derived by the 
fusion hierarhy (2^) and for e.g. 1 = 1,1 = 2 they are respectively 

l~l k-l 

hk{X) = — : n (2A + 2im) [] (2A + 2i + 2m) (2.15) 

2A + I n=0 

and 

l~l k-l 

hk{X) = Y[{2X + 2im) n (2A + 2i + 2m). (2.16) 

m=k n=0 

From the analyticity of the eigenvalues we obtain the Bethe ansatz equations 
e;+i(A„)e;i(A„)e,(A« - fi)2^e,(A« + l^)'^ei(A„) 

M 

= -ne2(A„-A^)e2(A, + A;3), (2.17) 

/3=1 

where = 2^"^ + 1. Recall that for / = 2s, E?^{Xq) is the permutation operator, 
therefore we obtain a local Hamiltonian 

;^log(t^^(A))|A=±f.+Ao (2.18) 
2ti aX 

and the corresponding energy eigenvalues 

1 2 M ^ 

^ ^ ~2^ S,?i (A, + (-)"^^ + f )(A, + (-)"fi - f ) • ^^-^^^ 



2.1.2 The boundary S matrix 

Now we can study the scattering of the excitations with the boundaries and determine 
the boundary S matrix. Before we do that it is necessary to determine what is the 
ground state and the low lying excitations for the model. The ground state, as in the 
^For general I 

l-l k-l 

hk{X) oc J]^ (2A + 2im) J| (2A + 2i + 2m). (2.14) 

m—k n—0 



8 



bulk, is the filled Dirac sea with strings of length q = 2s. The Bethe ansatz equations 
for the ground state become 

= n EggiX^ - \p)E,q{K + (2.20) 

/3=1 

Let us consider — and define 

^nm(A) = e|„_m+i|(A)e|„_m+3|(A) . . . e(„+m--3)('^)e(ji+m-l)(A) 

Km(A) = e|„_„|(A)ef„_^+2|('^)•••e(„+^_2)(A)e(„+„^)(A). (2.21) 

For the next we need the following notations 

X d 

an(A) = — — 2loge„(A) (2.22) 
and the Fourier transform of a„, is given by 

a„(a;) = e"'^. (2.23) 
We also need the following expressions 

(Znm{X). Anm(A)) = ^ log (X„^ ( A) , Km(A)) , (2.24) 

where their Fourier transforms are 

max(ri,m)a; i / / / \ \ i \ 

e 2 sinh (min(n,m))T7 

^nm(t^) = 2 coth(-)e 2 sinh(^(min(n, m))—j - 5nm- (2.26) 
Finally, the energy ( p.l9|) takes the form 



2 M 



E = -Y.llZqq{K + {-T^). (2.27) 



n=l Q=l 

The density for the ground state is given by 

ao(A) = Z,,{\ + n) + Z,,{\-n)~{A,,*a^){\) 

+ -^(-Z,,-(A)-Z,,+ (A) + a,(A) + A,,(A)) (2.28) 

where L = 2N is the length of the chain and * denotes the convolution. The solution 
of the last equation is 

ao(A) = e(A) + j{r{X,t) + r(A,r) + QW) (2.29) 
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where 

l + v4gg(a;) l + y4gq(a;) l + ylgg(a;) 

and 

<A)=E.(A-(-m, 5M = -^2|^. (2.31) 

We can write the above Fourier transforms in terms of trigonometric functions using 
the definitions of fig and A^g, 

- , , 6 2 —62 COShl^^l 1 . . 



2cosh(|)sinh(f ) ' > ' ' 2 coth(f ) sinh( 



and for > g — 1 



g-(x±-g)f 



S[UJ 



= E , w . ..^.x • (2.34) 



2cosh(f)' ^2cosh(7r(A-(-)*^]) 



Let us consider the state with u holes in the q sea, where u an even number. Then, 
in the thermodynamic hmit we obtain the density of the state form the Bethe ansatz 
equation, namely 

a(A) = (To(A) + ^X^(Ki(A-A,) + iri(A + A,)). (2.35) 
^ i=i 

The energy of the state with u holes in the q sea ( p. 27] ) is given by 

V 

E = Eo + J2<>^a), (2.36) 

0=1 

where Eo is the energy of the ground state and e(A) is the energy of the hole in the 
q sea. Finally, we compute the spin of the holes, and we can see from (|1.15|) that the 
spin of a hole in the q sea is ( [L.15 ) s = |, where / is an overall factor (see also [Q)- 
We consider the spin of the hole to be | for what follows. We conclude that the hole 
in the q sea is a particle like excitation with energy e, momentum p (e(A) = ~;^p(A)) 

= I 2cosh(.A = ^3 \| 5 '--'(-"(.A - (-)"n))(2,37) 

and spin s = | . We can easily check that in the scaling limit where X << Q the energy 
and momentum become from ( |2.37| ), 

e(A) ~ 2e-^^ cosh(7rA), p(A) ~ 2e""^ sinh(7rA) . (2.38) 
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Note that the excitations in the scahng hmit, as in the bulk case satisfy a massive 
relativistic dispersion relation. 

Having studied the excitations of the model we are ready to compute the complete 
boundary ^'-matrix. To do so we follow the formulation developed by Korepin, and 



later by Andrei and Destri , ||3^ . First we have to consider the so called quantization 
condition. 

(e2'^P5-l)|Ai,A2) = (2.39) 

where p is the momentum of the particle, the hole in our case. For the case of u (even) 
holes in the q sea, we compare the integrated density (|2.35|) with the quantization 
condition (|2.39| ). Having also in mind that, 

6(A) = --^ p{X) (2.40) 

we end up with the following expression for the boundary scattering amplitudes (the 
boundary 5" matrix will have the form = diag{a^, P^))^ 

a+a- = exp{2piL£' dX{(r{\) - e(A))} (2.41) 

and 

«^(A,?^) = fiX)ko{\)koiX)hiX,C^) (2.42) 

where 

A:o(Ai) = exp{2m£' \{YXKi{X - A,) + K,{X + A,)) + QiX))dX} (2.43) 
the X dependent part is 

h{Xi,^^) = exp{27rt j\{X,C^)dX}. (2.44) 

We are interested in the limit that s —>■ oo, in this limit we can easily see from ( |2.30| ) 
(see also |1^) that 

Qi.) KM, KA.) ^ H., e) = (2.45) 

with — g to be a fixed number as g oo and /(A) is an overall CDD factor given 
by 

/(A) =exp{ r e-^-47i7Y-} =tanhf(A + ^). (2.46) 
'^j-oo 2cosh(|-) uj ^ 2 2 

Using the above Fourier transforms we end up with the following expressions, 

, r f°° du ,,sinh(%^)e"^ ^ 

A;o A = exp{2 / — sinh 2za;A . . }, 
JO uj smh(2a;) 

h{X,e) = exp{-/ — sinh(2za;A)— — -}, (2.47) 

JO UJ cosh(ci;) J 
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where 2^^ = 2^^ — 2s + 1 is the renormahzed boundary parameter and A;o(A), ki{X,^^) 
are the ^ independent and the ^ dependent part, respectively of the usual XXX 
{SU{2)) reflection matrix (see e.g. [0). We notice from ( p.42| ) that there are two 
copies of the ^ independent part, whereas just one copy of the ^ dependent part exists. 
Recall, that we considered the diagonal K matrix ( |2.5|) in order to construct our chain, 
therefore we need to determine the other element of the boundary matrix. To do so 
we exploit the "duality" symmetry of the transfer matrix for —■6^- Hence, 

the other element is given by 

/3±(A,e^) = e2^-±_i(A)a=*=(A,|±), (2.48) 



the ratio ^,[^'-,] is the same as in the XXX model USSl but with a renormahzed 

a±(A,5±) " " 

boundary parameter. 



The above equations ( |2.42| )-( p.48|) , are simply combined to give the boundary 5" 



matrix as the tensor product of two rational {XXX) boundary S matrices (up to 
an overall CDD factor) i.e. 

S{X) = f{\)Ssui2){X,Ct)<^Ssui2){X,^t), (2.49) 
where = oo. This matrix coincides with the one found for the SU{2) 



PCM 1^ — in our case is a free parameter. 

However, in analogy with the bulk case, for finite s the boundary S matrix is 
expected to be of the form Ssui2){X,^'^) 'S> Srsos{s), where Srsos{s) is the bound- 
ary S matrix of the RSOS model. The Srsos matrix is a solution of the bound- 
ary Yang-Baxter equation in the RSOS representation (see e.g. [^-[^). The 



Ssu(2){X,C'^) 'S> Srsos{s) matrix should also describe the WZWk=2s boundary scat- 
tering and it should presumably reduce to the S matrix we found ( |2.49| ), as s — oo. 



2.2 The alternating s^, open spin chain 
2.2.1 The Bethe ansatz equations 

The corresponding transfer matrix t{X) for the open alternating chain of 2N sites and 
= y, = ^ spins {q^ — q'^ = 1) is (see also e.g., ||2^, [§]), 



t(A) = tro K^{\) To{X) K^iX) fo(A) , (2.50) 
where To(A) is the monodromy matrix defined previously in ( [1.21| ) and 



fo(A) = Rl,{X - Q)Rlo{X + Q)--- /?^;v-io(A - ^)RlNoiX + ^) ■ (2-51) 
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Again here we use the diagonal K matrix ( |2.5| ) and the fusion hierarchy ( |2.6| ) with 

5[T{\)] = Cl(A-^] + ^)^C2(A + ^^ + ^)'^, 

5[f{X)] = Ci(A + fi + 2)^C2(A-^^ + «)'^ (2.52) 

and 

C„(A) = (-A + is'' + ^)(A + is'' + '-). (2.53) 

The eigenvalues of the transfer matrix, after we fuse the auxiliary space, are given by 



Ap'')(A) = 5:/i,(A)^(A)a«(A-^])^4^^(A + ^])^af (A-fi)^af (A + fi) 

(A - A,- - f ) (A - A,- + z/ + i) 



N 



k=0 

M 

'1 2' 



n (A - A,- + i{k - 1) + f ) (A - A, +ik + ^) 
n + (A + A,+z/ + f) 

/ii (A + A,- + i{k - 1) + 1) (A + A,- + zfc + I) ^ ' ' 

provided that {A-,} satisfy the Bethe ansatz equations 

n en{\a - f^) V(Aa + f])^ei(A,)e;+^ (A„)e;_^ (A„) 
= - n e2(A, - A^)e2(A, + A^). (2.55) 

/3=1 

For I = 2s\ R^{Xq) becomes the permutation operator, therefore we can obtain a local 
Hamiltonian for the open chain 

H oc ^|^log(^'(A))lA=o+Ai + £:|^l°g(^'(^))lA=-n+Ag (2.56) 



and 



1 q'' 



EE 



4^nt^i,tl (A,-fi + ^)(A,-f]-^) 

~ ^n?i§ (A, + f^+¥)(A, + f^-¥)' ^^'^^^ 

2.2.2 The boundary 5 matrix 

The ground state consists of two filled Dirac seas with strings of length = 2s"', 
n = 1,2. Then the Bethe ansatz equations for the ground state become 



n n EnAK - K)EnAK + A^,) (2.58) 
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where n can be q^, and let us consider for simplicity > g^. Finally, the energy 
( p.571 ) takes the form 



q2 M, ^ g2 M, 



^ = E E + ^) - ^ E E - (2.59) 

The density that describes the ground state is given by the following integral equations, 

2 2 

<W = E U^n,{X + Z„,(A + n))-Y: (^n, * a^)(A) 

-1 i6 .1 

1 

+ -(-Z„,+ (A) - Z„,-(A) + a„(A) + J2 ^nM)) (2.60) 

where n can be or g^. The solution of the above integral equation is given by 

ao"(A) = e"(A) + i(r"(A, t) + r"(A, T) + Q"(A)) (2.61) 

where 

9^ g2 



f"(a;,e^) = -^Z„±Hi?„,H, Kr(^) = E ^™(^)^™(^) (2-62) 



and 



e"(A) = ^E^(A-(-)^l]), (2.63) 
^ i=i 

and its Fourier transform is 

S"(a;) = 5: (ZMico). (2.64) 
Here i? is the inverse of the kernel K of the system of the linear equations (|2.6CI| ), 

Knm{i^) = (1 + Anm{t^))6nm + -4„m(u;)(l - 5„„) (2.65) 
1 . 

Rnmiyj) = ^ E ((1 + Ajj{uj))6nm{l - ^nj) - Anm(^)(l " ^nm)), (2.66) 

cletK j^gi 

where the determinant of K is, in terms of trigonometric functions, 

detK = 4coth2(|)e~^ sinh(g2|) sinh(|). (2.67) 

In particular, K^"^ has the following explicit form in terms of trigonometric functions 

e-t ^y,,,, sinh((g^-2)f 



2cosh(f) 2cosh(f)sinh((gi-l)f; 
= Af(-0 = ^^ (2.68) 

14 



gi(a;) = l±^, Q\uj) = k^\uj) + \-- + '^""^^^^ , . (2.69) 

^ ' 2cosh(f)' ^ ' ' 2cosh(f) 2cosh(f)sinh(2^) ^ ' 



The ^ dependent part for > is 



finally, 



^"(^) = E— — TT^— T^- (2-71 



2cosh(f) tiAcos\i[Ti{\ + {-yn) 



Let us consider the state with i/„ holes in the g" sea, where is an even number. 
Then in the thermodynamic limit we obtain the density of the state from the Bethe 
ansatz equations, namely 



1 

a"(A) = <(A) + y E(^r(^ - ^i) + ^r(^ + ^^)) (2-72) 
^ i=i 

the energy of the state with Vn holes in the g" seas is given ( |2.59| ) by 

E = Eo + Y.^^iK). (2.73) 

a=l 

where Eq is the energy of the ground state and e"(A) is the energy of the hole in the 
sea. Finally, we compute the spin of the holes from ( |1.28|) , and we can verify that 

the spin of a hole in the sea is = | whereas the spin of a hole in the q^ sea is 
= 0. We conclude that the hole in the g" sea is a particle like excitation with energy 

e", momentum (e"(A) = ^^p"(A)) 

e"(A) = E 



1=1 



4cosh(7r(A- (-)'!]) 



P"(A) = ±^ + E^tan-l(sinh7^(A-(-)'^])), (2.74) 
^ 1=1 ^ 

and spin = |, = 0. We can easily check that in the scaling limit, A << Q the 
energy and momentum become from ( p.74| ), 

e'^(A) ~ e""^ cosh(7rA) , p"(A) ~ e"^^ sinh(7rA) , (2.75) 

the factor e~'"^ provides a mass scale for the system. Note that in the presence of 
boundaries the excitations, in the scaling limit, satisfy a massive relativistic dispersion 
relation (|2.75| ) whereas in the bulk case the excitations are massless relativistic particles 
( |1.29| ) . This is a very interesting phenomenon which is presumable related to the type 



of the boundaries. The boundaries we impose, force a left (right) mover to reflect as 
a left (right) mover. It is possible that the boundary can reflect a left mover to a 
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right one and vice versa. This type of boundaries would probably lead to massless 
excitations in the scaling limit. 

Again, we consider the quantization condition in order to compute the 

exact reflection matrices, namely 

(e'*^^"^?"- l)|Ai,A2) = (2.76) 

where is the momentum of the hole in the sea. For the case of z^" holes in the 
sea, we compare the integrated density (|2.72| ) with the quantization condition ( |2.76| ). 
Having also in mind that, 

e"(A) = p-{X) (2.77) 
we end up with the following expression for the boundary scattering amplitudes (S'^ = 



exp[2piL j^' d\[a{\) - e(A))} (2.78) 



and 

al{\e) = k-,{X)k-,{\e) (2.79) 

where 

A;^(Ai) = exp{7r2 ^(ir""(A - A^) + K""(A + A^)) + Q"(A)(iA} (2.80) 
•^0 i=l 

the X dependent part is 

A;^(Ai, e^) = exp{27r? j^' r"(A, £,^)d\]. (2.81) 
We are interested in the limit that g" — oo, in this limit we can easily verify that 

and 

with — to be a fixed number as g" — >■ oo. Using the above Fourier transforms we 
end up with the following expressions, 

f du , , sinh(%^)e"^ 1 

kUX,0 = exp{-/ _sinh(2z..A)— — -}, (2.84) 
^ JO uj cosh(a;) 
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where 2it = 2^^ -25^ + 1, 2|^ oo, are the renomarhzed boundary parameters, 
fcQ(A) and /c"(A) are the ^ independent and the ^ dependent part, respectively of the 
usual XXX {SU{2)) reflection matrix (see e.g. |38|). Exactly as in the case of the 
spin s open chain, we consider the diagonal K matrix, therefore we need to determine 
the other element of each boundary matrix. We exploit the "duality" symmetry 
of the transfer matrix for and we find that the other diagonal element is 

given by 

m\^n) = e2.-±-i(AX(A,|^). (2.85) 



The ratio ^f.,^^" is the same as in the XXX model but with a renormalized boundary 
parameter. We observe for the alternating chain as well as in the spin s chain that 
only one free boundary parameter (^^ oo) survives. 

Two copies of the rational {XXX) reflection matrix were computed, one for each 
excitation. Therefore, we conclude that the boundary S matrix should be of the 
structure, 

SiX) = Ssu(2)iX,it)^Ssui2)iX,^t)- (2.86) 

We assume that this matrix should also coincide with the one of the PC Mi. A cal- 
culation of the boundary S matrix from the field theory point of view would probably 
confirm our results. We have to mention that there have been some studies for the 
SU{2) PCM with WZW term with boundaries, but mainly in the context of 



quantum impurity (Kondo) problem. In particular, in |Q the authors considered dy- 
namical boundaries, i.e. they considered quantum impurities at the boundaries, and 
they derived the corresponding reflection matrices. 



3 Discussion 

The XXX spin s and the alternating s^, (s^ ~ ^'^ — \) chains were explored. 
For both models the Bethe ansatz equations were derived using fusion, and the ex- 
act boundary 5" matrices were computed. We were particularly interested in the case 
that s, oo. More speciflcally, for the spin s — oo chain the boundary scat- 

tering amplitudes were simply combined to give the boundary S matrix of the form 

/(A)55C/(2)(A,|±)055C/{2)(A,|2'=). 

For the alternating spin chain two different types of excitations exist: \ and 
spin respectively. The boundary scattering for each excitation was studied and the 
corresponding reflection matrices were derived. Two copies of the XXX boundary S 
matrix were computed (s* — >■ oo), and the boundary S matrix was given as a tensor 
product of two rational boundary S matrices. This matrix is also expected to coincide 
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with the one of the SU(2) PC Mi. Note that we could end up with to same result if 
we started from the anisotropic spin chains and then take the isotropic and s, ^ oo 
limits. 

It would be also interesting to consider dynamical K matrices in order to 
construct the open spin chain and then to study the reflection of the particle-like 
excitations with the dynamical boundary. Another interesting aspect would be the 
study of the thermodynamics of the alternating spin chain via the TBA. The main 
purpose would be the derivation of the central charge (see e.g. |2^, |2^), of the model 
for both bulk and boundary cases. We hope to report on these issues in a future 
pubhcation. 
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